%%
load c:\manoj\projects\ace\pred_julia ace_matrix;
load c:\manoj\projects\ace\Julia_W_new;
i=find(floor(fday(:,1))==datenum(datenum(2006,10,01)-datenum(2000,1,1)));
subplot('position',[0.15,0.15,0.8,0.2]);
plot(time_ax+datenum(2000,1,1),ace_matrix(i,Julia_W(i).k),'k','LineWidth',2);
set(gca,'FontSize',16);
datetick;
axis([datenum(2000,1,1)+Julia_W(i).fday+(8+5)/24 datenum(2000,1,1)+Julia_W(i).fday+(16+5)/24 -inf inf]);
grid on;
xlabel('UT (hours)');
ylabel('IEF Ey mV/m');


subplot('position',[0.15,0.4,0.8,0.5]);
plot(time_ax+datenum(2000,1,1),w(i,Julia_W(i).k),'b','LineWidth',2);
set(gca,'FontSize',16);
hold on;
plot(time_ax+datenum(2000,1,1),julia_model_drift(i,Julia_W(i).k),'r','LineWidth',2);
plot(time_ax+datenum(2000,1,1),w(i,Julia_W(i).k)-julia_model_drift(i,Julia_W(i).k),'b-');
datetick;
axis([datenum(2000,1,1)+Julia_W(i).fday+(8+5)/24 datenum(2000,1,1)+Julia_W(i).fday+(16+5)/24 -inf inf]);
grid on;
set(gca,'XTickLabel',[]);
date_string = datestr(time_ax(1)+datenum(2000,1,1),29);
title(sprintf('Date - %s',date_string));
ylabel('V_z (m/s)');


